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Abstract 



We consider a system of differential equations the behavior of which 
solutions possesses several properties characteristic of the blood pressure 
distribution. The system can be used for a compartmental modeling [8lll5j 
I of the cardiovascular system. It admits a unique bounded solution such 

that all coordinates of the solution are separated from zero by positive 
numbers, and which is periodic, eventually periodic or almost periodic 
^ , depending on the moments of heart contraction. Appropriate numerical 

H ' simulations are provided. 



1 Introduction and Preliminaries 

A living organism can be considered as a complex system of nonlinear oscillating 
structures of different origins. The system as an integrity is a "constellation" 
of all of its oscillations. The oscillators are connected with each other to build 
oscillating chains. The essential task of controlling processes in a living organism 
is to sustain the oscillating activity of certain significantly unstable elements in 
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such a way, that the amphtudes of their oscillations don't become excessively 
large or excessively small. It is, therefore, natural to assume that the data 
of oscillations should be considered as one of the most important sources of 
information about the condition of the organism. 

The movement of any fluid in a certain direction is the result of pressure 
differences in the starting and ending points of the system. The difference of 
pressure in heart chambers and in aorta, in aorta and outgoing blood vessels is 
the cause of blood flow in an organism. Last decades there has been a rise of 
intensive investigations of blood pressure in its connections with other parameters 
of the cardiovascular system [H [3],[7j-[n],[lll [HI [2Q],[23]-[25],[29]. At the same 
time we should recognize that due to the nonlinear properties of blood vessel 
walls and their variability along the length of a vessel the relation between actual 
oscillations of pressure and of blood volume [10] , [20] , [25] is very complex, and 
consequently, it is reasonable to consider the distribution of blood pressure out 
of it's connection with other parameters, particularly blood volume. There are 
many different approaches to the problem of modeling cardiovascular dynamics 
including, for example, modeling of the coupling mechanisms as a system of 
differential equations with delay [HI [I2] or as a system of ordinary differential 
equations [TH UHl [29] . Compartmental models have been effectively applied to 
investigate different type processes in chemistry [22], medicine [HI [121 [13 [16] . 
epidemiology [21], ecology [19], pharmacokinetics [6l [8], [Ml [28] • The idea of 
compartments makes our investigation more abstract. Particularly, we do not 
specify the physiological meaning of the compartments nor do we specify the 
number of these parts. Only compartments which have an essential influence on 
the aortic pressure are considered. 

It is obvious that all parts of the cardiovascular system should be involved if, 
let say, the blood flow problem is investigated, because the total volume of the 
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blood is constant in the time, but the pressure could be discussed locally, as, for 
example, in [III [18], where blood pressure for only compartment is studied. In 
our paper we focus on the interactions of the blood pressures of several parts of 
the system, but not all of them. We suppose that there exists, as a consequence 
of the blood float, a mutual interaction of blood pressures in different parts of 
the system. One of them, the systemic arterial pressure, should be singled out, 
as it is disturbed impulsively with large simultaneous inflow of blood in the time 
of heart contraction. Other compartments of the system do not directly undergo 
the impulsive action, so their change is continuous. We want to note that in this 
paper we consider exclusively the aorta and the parts of the cardiovascular system 
interacting with the aorta through blood pressure, that is, the compartments 
belonging to the systemic circulation only are considered. The blood pressure in 
pulmonary artery also changes impulsively, and the differential equations system 
we propose can therefore be also used to model the pressure in this artery and 
its neighboring compartments. We, however, leave it out of our discussion. 

The goal of our paper is not to build a final model of cardiovascular system, 
but rather to introduce a system of differential equations, in which one coordinate 
is perturbed impulsively and other coordinates exhibit a continuous change, thus 
having certain properties characteristic of the actual blood pressure behavior. 
The obtained results now need qualitative experimental confirmation. We hope 
that our results will serve as a basis for further experimental and clinical studies, 
which can lead to better understanding of the cardiovascular system. 

Let Z, N and M be the sets of all integers, natural and real numbers, respec- 
tively, ]R+ = [0,00), and let || ■ || be the Euclidean norm in M", n eN. 

We denote Cq the aorta and arteries, and Pq blood pressure in Cq. Let Ci,i = 
l,m, be other vessels and organs, which adjoin to Co, and have an essential 
influence on fluctuations of Pq, and Pi,i = l,m, are the blood pressure values in 
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compartments Cj, 1, m. 

In our paper we suppose that the pressure variables satisfy the following 
system of differential equations with impulses 

dPoit) 



dt 
dPijt) 

dt 
dP2{t) 

dt 



-koPo — go{Po — Pi, Pq — P2, ■ ■ ■ ,Po — Pm), 
-hP, + g,{PQ-P^), 

-k2P2 + g2iPo-P2), 



^Prn{t) _ 1 p , (p _ p \ 

APo\t=e, = Io + MPo), (1) 

where APo|t=e, = Po{Oi+) - Po{9i) , Po{9i+) = \imt^e,+ Po{t) . 

The following assumptions for ([T]) will be needed throughout the paper: 



(CI) real constants Iq, ki, i = 0,m, are positive; 



(C2) Jq, gi are real valued functions, Jo(0) = 0,gi{0) = 0,i = l,m,Jo{z) > 
0, gi{z) > 0, if 2; > 0, go{0, 0, . . . , 0) = 0, 5^0(^1, ^2, • • • , ^m) > 0, if Zj > 



0,j = l,m; 



(C3) the functions Jo,gi,0,m, satisfy the Lipschitz condition 



i{z^) - gi{z'^)\ < li\z^ - z^\,i = l,m, 



\go{zi, . . . , Zm) - go{wi,...,Wm)\ <^ojYl 



\Zi — w. 



1=1, m 



\J,(z')-J,{z')\<lj\z'-z'\- (2) 
(C4) there exists a number u > 0, such that iu < 6i < {i + l)uj, i E Z; 



(C5) there exist positive real constants mj, z = 0, m, rrij such that 
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sup gi{z) =mi,i = l,m, sup go{zi, . . . , z^) = nio 

z>0 Zj>0,j=l,m 

sup Jo{z) = mj. (3) 

2>0 

We assume that atmospheric pressure has the zero value. Since pressures can 
not be negative, the main goal of our paper is to find the conditions which guar- 
antee the existence of bounded positive solutions. Moreover, the coordinates of 
the solutions are separated from zero by some positive constants, and eventually 
the first coordinate must be larger than any another. The solutions should be 
exponentially stable under certain conditions, and their various oscillatory prop- 
erties (periodicity, eventually periodicity, almost periodicity) correspondly to the 
particular properties of the sequence of moments 6'j. 

Differential equations ([H) as a model of the physiological process are developed 
using models, which can be found, for example, in [H [9l [HI [HI [2l] and many 
others. To clarify that, we shall describe the two following models. 

In paper [18] the following differential equation is considered for the peripheric 
blood pressure Pp 

where C is the total arterial compliance, R— the arterial resistence, Q{t)— the 
continuous blood flow. The formula for the continuous blood flow evaluation is 
given by Q{t) = E°M}:ifpMl^ where Pq— is the systemic arterial pressure, and r— 
is the aortic impedence. As we suppose that the systemic arterial pressure is 
larger than the blood periferic pressure, Q{t) is a positive function. 

For the systemic arterial pressure the following differential equation with 
impulses has been discussed in book [H], 

dPo 1 



dt RC 



P,{e,+)-Po{9,) = ^, (5) 



where 9i = iT,i G Z, are prescribed moments, Po(^i+) is the right hmit's value, 
Vj are the stroke volumes, and C is the systemic arteria compliance. 

It is natural to consider cardiovascular problems using discontinuous dynam- 
ics theory. One of the interesting approaches which should be mentioned is the 
method of circle mappings [71 [TT] for cardiac arrhythmias. In [3] we consider the 
extended form of equation ([5]) for the systemic arterial pressure, and we investi- 
gate the problem of existence of positive periodic and almost periodic solutions, 
oscillations of the equation. The case when the moments of heart contraction are 
not prescribed, but are caused by a certain value of the systemic arterial pressure 
has also been investigated. 

We have further developed the models of our predecessors in several ways: 

1) The nonlinear perturbations caused by interactions are used in the model, 
both in differential equations and the impulsive part of the model. Function 
9o{Po — Pi, Pq — P2, ■ ■ ■ , Pq — Pm) is positive since the outflow from the aorta to 
neighbouring compartments reduces the systemic arterial pressure, all differences 
in the function are assumed to be positive. Nonlinearities gi{PQ — Pi),i = 

are positive, when the differences again are positive, since the differences generate 
inflow into compartments, which naturally implies the rise of pressure. Every 
function gi{Po — Pi) does not depend on Pj,j 7^ i, since we assume that there is 
no essential mutual interaction between compatments; 

2) The equations for the compartments are derived from equation (jl]) for 
peripheric blood pressure. The nonnegative differences ^"^--P^W g^^.^ replaced by 
the nonlinear functions gi{Po — Pi); 

3) Jo is a constant part of the instanteneous change of systemic arterial pres- 
sure; 

4) Condition (C5) is chosen, assuming that the results obtained in our investi- 
gation can be interpreted for the blood pressure system activated by an artificial 



6 



pacemaker, or for the case when fluctuations are regular. That is, we exclude 
the case when the set of points of discontinuity may have so called accumulation 
points. 

We must also note that our results can be easily reconsidered when condition 
(C5) is replaced by condition < /i < 6'j+i — < /ii < oo or by the condition of 
existence of lim sup j.^^^ <*o,^+r) _ -y^^j^gi-g ^(to, h+T) is the number of elements 
0i in the interval (to, ^0 + ^)1 and p is a nonnegative real number, uniform for all 
to G M. 

The mathematical background of our investigation, the theory of impulsive 
differential equations, has been developed in [2], HI [5], [131 [13 E?]. The results of 
the paper concerning eventual periodicity, positiveness of solutions, domination 
of a coordinate over all other coordinates, are new. The role of the impulsive 
action as the one determining the relations is investigated for the first time. We 
are aware that in reality the role of the impulses is diminished, for example, 
by delays, and we plan to continue the investigation of the interactions of blood 
pressure levels in different parts of the cardiovascular system in our future papers. 

The paper is organized in the following manner. In the next section we find 
the conditions for the equation such that there exists a unique solution bounded 
on the whole M . The solution is periodic, almost periodic, or eventually periodic 
if the sequence of discontinuity moments has an appropriate property. Section 
3 contains the results on the stability and positiveness of the solutions. In the 
last section we specify the choice of the moments of discontinuity such that they 
may be determined artificially by a particular map, and by the initial moment. 
In this case the approach can be considered as a way to discuss the pacemaker's 
design. Appropriate numerical simulations are provided. 
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2 Bounded, periodic, eventually periodic and 
almost periodic solutions 

In this section we shall prove that the system admits oscillating solutions. The 
existence of an eventually periodic solution is proved for the first time. The as- 
sertions about the discontinuous periodic and almost periodic solutions are spec- 
ified versions of general theorems of the theory of impulsive differential equations 

Firstly, let us define the solutions of system ([1]). 

A left continuous function P{t) : R R^^^ is in the set VC{R), if it is 
continuous on R, except at the points 6i,i G Z, where its first coordinate may 
have discontinuities of the first kind, and is left continuous. 

A function P{t) e VC{R) is a solution of (P) if: 

(1) the differential equation is satisfied for P{t),t G R, except at the points 
6i,i & Z, where it holds for the left derivative of P{t); 

(2) the jump equation is satisfied by Po(t) for every i E Z. 
We shall also need the following set of functions. 

A left continuous function P(t) is in the set of functions VC{R+),Rj^ = 
[to, oo), to £ R, if it is continuous on [to, oo) except at the points 9i, < 6i < oo, 
where its first coordinate may have discontinuities of the first kind. 

A solution P{t) = {Po(t), Pi{t), . . . , Pm(f)) of ([T]) on [to, oo) is a function in 
VC{R+) such that: 

(1) the differential equation is satisfied for P{t) on except at the points 
6i,9i > to, where it holds for the left derivative of P{t); 

(3) the jump equation is satisfied by Po(t) for every i, 9i > to- 
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One can easily show, using the theory of impulsive differential equations 
[ITlETj, that under the aforementioned conditions a solution P{t, to, ttq), (to, vq) G 
R X R^'^^, ttq = (tTq, tTq, . . . , TT^), of ([1]) exists and is unique on [to, oo). 

One can also find that the solution satisfies the following integral equation 

t 

j e"^°(*-^)(?o(Po(s) - Pi{s),Po{s) - P2{s), Po{s) - PUs))ds + 

to 

to<ei<t 

P,{t) = e-^'(*-*«)7r^ + 
t 

J e-'=^(*-^)(/,(Fo(s) - P^{s))ds, I = M^. (6) 

to 

We say that a sequence 6i has the p— property, p E N,ii Oi^p = Oi+puj for all 
integers i. Moreover, we say that the sequence has the p— property eventually, if 
Oi^p — 9i ^ pu as i ^ oo. 

Consider a strictly ordered sequence of real numbers ti,i G Z. Denote = 
ti+j — ti,i,j G Z, and define the sequences {fj}i,i,j G Z. Following [271 [H] we 
call this family of sequences equipotentially almost periodic if for an arbitrary 
positive e there exists a relatively dense set of e— almost periods, common for all 
sequences {t^j},j G Z. 

It is proved in |4| (see also [27]) that the family 6j, i,j G Z, is equipotentially 
almost periodic if the sequence 9i — iuj,i G Z, is almost periodic. That is, if we 
chose an almost periodic sequence G Z, such that < < u,i E Z, then 
the sequence 9i, where 6i = iu + E Z, is equipotentially almost periodic. 

Let [a,b] be an oriented interval, that is [a,b] = [a,b], if a < 6, and [a,b] = 
[b, a], a a > b. 

We say that a function (f){t) in VC{R) is eventually pcij— periodic if the se- 
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quence of discontinuities 9i has the p— property eventually, and 



t— »oo 



\im[(l){t + pu) - (j){t)] = 0, 



(7) 



for all t such that t ^ [9i, 9i+p — pu] fl [6'j+i, — pu]. 

Next, we shall give a definition from [271 [T3] which is slightly modified for 
our case. 

A function (p{t) in VC{R) is a discontinuous almost periodic (d.a.p.) function 



(a) the function is uniformly continuous on the union of all intervals of the 
continuity {9i,9i+i),i G Z; 

(b) for every positive e there exists a relatively dense set F of almost periods 
such that if 7 G r, then \\(j)(t + 7) — (pit)\\ < e for all t G -R such that 
\t - 9i\ > e,i e Z. 

We may assume that 



Theorem 2.1 If conditions (CI) — (C6) are fulfilled, then there exists a unique 
solution of (Cy bounded on R and: 

1. the hounded solution has the period puj if the sequence 9i has the p— property 
for a fixed p & N; 

2. the bounded solution is an eventually pu— periodic function if the sequence 
9i has the p— property eventually for a fixed p E N; 

3. the bounded solution is a d.a.p. function if the sequence 9i — iu,i G Z, is 
almost periodic. 



if: 
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Proof: We shall use the norm |(/)|o = snpp, \(f){t)\ for scalar valued functions 
defined on R, and ||0||o = sup^ ||(^(t)|| for vector-valued functions. Denote by 
Q a subset of VC{R) such that if G fi, = (0O) 0i) • • • 5 ^m), then |0o|o < 



(^ + T^), \c^.\o<f,^ = l,m. 

One can easily verify that ([T]) has a bounded solution P{t) if and only if P{t) 

is a bounded on R solution of the following integral equation 

t 

Po{t) = - j e-'''^^'~'^go{Po{s)-P,{s),...,Po{s)-PUs))ds + 

—00 

ei<t 

t 



P,{t) = j e-'='(*-^)(7,(Po(s) - P^{s))ds, i = l,m. (8) 

—00 

Define on Q an operator T such that 

t 

—00 

^e-'^°(*-^')(/o + Jo(0o(^.))), 

0i<t 

t 

{TcPUt)= j e-'^'^'-''>g,iMs)-Ms))ds,z = T-i^. (9) 

—00 

We can verify that T : Q ^ Q. Take (f),ip G Q. One can obtain that 

\\T(j)-T4j\\o<L{lo,h,...,lm,lj)\\<p-i^\\o, 

and, consequently, condition (C6) implies that the operator is contractive. It is 
not difficult to check that the space Q is complete. Thus, there exists a unique 
solution of ([1]) from Q. 

Assume now that the sequence 6i has the property. Then, using the stan- 
dard method, starting with a pu;— periodic function from Q, we can, using the 
operator T, construct a sequence of pcu— periodic approximations of the bounded 
solution, which is also pcu— periodic. 
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Consider the case where 6i has the property eventually. In order to prove 
assertion 2) of the theorem, it is sufficient, according to the above discussion of 
boundedness of the solution, to prove that Tcj) is eventually pcu— periodic, if (f) is 
an eventually pa;— periodic function. 

Fix a positive number e. Since 6i has the property eventually, and <p is 
eventually pcu— periodic function, one can find positive numbers T2 > Ti such 
that: 

1. if t > T2 and 9i + e < t < ftj+i — e for some i E Z, then 9i^p < t+puj < ^^j+i+p 
and \\(j){t + pLo) - 0(t)|| < e; 

2. 2e-'=o(^2-Ti)r!Bi + IIMS^] < ^ 2e-'^°(^2-^^)^ < e. 

L fco 1— C 0" J ' fci 

We have that 

rt+puj 

\{T<f))o{t+pu;)-{T(P)o{t)\ = I- / e-''<^^'+P^-'^go{Ms)-Ms), . . . ,Ms)-Ms))ds+ 



j2 e-^o(*+P--^^.)(/Q+j^(0^(^.)))+ r e-'^'^'-'^go{Ms)-Ms),...Ao{s)-<l>m{s))ds- 

^j- 1 . —00 



rTi 



e ^«(*"^)2mods+ 



^ e-^o(*-^')2mj + f e-'''^'-'^lov^\\<j){s + puj) -(P{s)\\ds+ 



Ti<ei<t 0i<f 



Similarly, one can obtain that 

|(T0),(t + puj) - (T0),(t)| < e[l + -i + 4m,e^o^- ^ = 1, m. 

ki 1 — e"'''^ 

That is, there exists a unique eventually pcj— periodic solution of ([T]). 
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Now assume that the sequence 9i — iu is almost periodic. From the previous 
discussion we see that to prove the existence of a d.a.p. solution we need only 
to verify that the function T(p is d.a.p. if is a d.a.p. function. By Lemma 35 
[27j (see also [1]), for a given positive e there exist a real number u,0 < u < e, 
and relatively dense sets of real numbers F and integers H, such that: 

2. \9^-j\<u,keZ,heH,-fe F. 

Take the real number u such that — 0(^2)|| < e if ^1,^2 belong to the same 

interval of continuity of the function (f){t) and |ti — ^2! < J^- We have that if 
|t — > e,i & Z, then 

-00 

J2 e-^"(*+^-^»)(/o+Jo(0(^.)))- r e-'^°(*-^)<7o(0o(s)-0i(s), • • • , Ms)-Ms))ds- 
^e-'=°(*-^»)(/o + Jo(m)))l< r e-'=«(*-)/oV^||0(s + 7)-0(^)IM^+ 

J2^~'°^'~''^ij\\m+h) - m)\\ + E f'^^ e-'^^^'-^hmods < 



i<t ei<t ■ 



and 



Jo 2e'=»'' ^ . 



|(ri>).(t + 7) - < e|A + 4mie'--^^],i = l,m. 



That is, T(p is a d.a.p. function. The theorem is proved. 
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3 Stability and positiveness 

In this section we show that the bounded solution is exponentially stable, and 
that every coordinate of the solution is separated from zero by some positive 
number. Moreover, the first coordinate of the solution, which is the value of the 
systemic arterial pressure, is higher than any other pressure value. 

Let us denote the solution bounded on R as ^(t) = {Co{t) ^ ^lit) i ■ ■ ■ iCm(t))- 
Next, we find the conditions for the positiveness of this solution. Assume addi- 
tionally that 



(C8) ^ + ^<Io^,^ = l,m. 
Using (IHl) we obtain that 



and 



Hence, 



<T^,^ei?,z = l,m. 

rvi 



Ut) - m > loY^^ --^--^ = 5>0,teR,i = l,m, (10) 

and using again one can see that 

9i 



where gi is the minimal value of the function gi{z) for 2 G [S,^ + ^"l-koi^ + 
The following assertion is proved. 

Theorem 3.1 Assume that conditions (CI) — (C8) are valid. Then for the 



bounded solution ^{t) of ^ there exist positive constants = l,m, fii,i 



0, m + 1 such that the inequalities ^oit) — ^i{t) > Uj,^i(t) > yUj are valid for all 

1, j, and t E R. 
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Let us give the definition of uniform exponential stability of the solution. 
Denote P{t) = P(f,to, Po), a solution of (□). 

Definition 3.1 The solution ^{t) is called uniformly exponentially stable if 
there exists a number a ^ R, a > 0, such that for every e > there exists a num- 
ber 6 = 6{e) such that the inequality \\P(t) — < e exp(— a(t — to))i Vt > tQ, 
holds, z/||Po-^(to)|| <5. 

Fix a positive number fx, < o" < minj^o;^A;j, denote m{lQ,li, . . . J^Jj) = 
1 — max{^|^3^ + Y^^qi^p^) F^' ^ = 0; ""^l) and assume that the Lipschitz coef- 
ficients are small so that 

(C9) m{lo,h,...,l^,lj)>0. 

Theorem 3.2 Assume that conditions (CI) - (C7), (C9) are valid. Then the 
bounded solution ^{t) of (Qp is uniformly exponentially stable. 

Proof. One can see that v(t) = P{t) — ^{t),v = {vq, t>i, . . . , Vm), is a solution 
of the equation 

^ = -kovo + wo{v), 
at 

dv. 



— = -kiVi + Wi{vo,Vi),i = l,m,t^ Oj, 
Afo|t=ej = uo{vo). (11) 



where functions Wi,i = 0,m, and Uq satisfy the following conditions: 

\wo{v) I < ZoV^2m| It^l I; 
\wi{vi,Vo)\ < li{\vi\ + \vo\),i G Z; 
\uo{vo)\ < lj\vo\. 

Thus, the problem of the stability of ^o(^) is reduced to the stability of the zero 
solution of (ITTl) . Fix e > and denote K = KUq, h, . . . , Im, h, = i, , ^ , — tt; 
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where 6 & R, 6 > 0, and take 6 so small that K{lo, li, . . . ,lrn,h,S) < e- Assume, 
without loss of generality, that to = 0. Let v{t, vq) be a solution of ( fTTi) such that 
v{0,Vo) = vo= {v^,v\ . . .^v""). 

Denote by a set of all functions i/j = {i/jq, ipi,'ilj2, ■ ■ ■ , V'm), defined on = 
[0,oo), such that: 1)^(0) = Vo] 2) e VC{R+)] 3) < K exp{-at) if 



t > 0,i = 0,m. Define an operator 11 on \1/ such that if ■?/' G then 



t 



(n^)o = e-'^oV + fe-'"^'-'^woWs))ds+ e'''^'~''^Mvo{0^))), 
o<ei<t 
t 

(nV). = e-^»V+ [ e~''^^'-'^w,{vo{s),v,{s))ds,t = l~^. (12) 







We shall show that 11 : ^' — > Indeed, for t > it is true that 

\{IIiIj)o\ < exp{—kot)6 + / exp{—ko{t — s))loV2mKK exp{—as)ds+ 



1 \/'?m 1 p(^0—cr)ul 

o<ei<t ° 



Similarly, 



|(nV')i| < exp{-at)[5 + 2K ^ ] < Kexp{-at),i = l,m. 
Let ipi,ip2 € ^r?- Then, we have that 

I in?/'! -nV'2|| < iv(/o,/l,...,/m,/j)sUp||V'l -^2||, 

i>0 

where the coefficient is described in (C7). Using the contraction mapping ar- 
gument one can conclude that there exists a unique fixed point v{t, Vq) of the 
operator 11 : \1/ — \I' which is a solution of ( ITT]) . Theorem is proved. 

From ( fTOl) it follows that the bounded solution has the first coordinate larger 
than any other coordinate, and the attractiveness of the solution implies that 
any other solution in its neighbourhood eventually has the first coordinate as 
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its largest coordinate. Now, we can make some physiological conclusions: the 
normal state of the distribution of blood pressure is such that systemic arterial 
pressure is higher than any other pressure. The initial state may be odd, that 
is, the initial systemic arterial pressure may be lower than the pressure in a 
neighbouring compartment, but after a certain period of time the state becomes 
normal. We may call this period the time of stabilization of a solution. 

In order to carry out numerical simulations of the obtained theoretical results 
we consider the following equation 

F^ = -0.5Fo-0.1(Po-^'i), 
P[ = -0,7P, + 0.1{Po-Pi), 
P^ = -1.2Pi + 0.1{Po-P2), 

APo|t=e. = 0.07, (13) 

where 6i = i + sin(i) —sm(\^i)\. It is proved in [4J that the sequence ^| sin(i) — 
sm(\^i)\ is almost periodic. The theoretical part of the paper implies that 
system (1131) has a discontinuous almost periodic solution. 

As can be seen from Fig. 1, the solution P{t),P{0) = (0.08,0.05,0.03) is 
approaching the almost periodic solution as time increases. 



17 



0.18 






0.16 






0.14 






0.12 




III 


0.1 






0.08 






0.06 






0.04 




0.02 







1 1 1 1 1 1 



5 10 15 20 25 30 35 

time t 



Figure 1: The graph of the first coordinate Po{t) of the solution is shown in red, 
of the second coordinate Pi{t) in blue, and the third coordinate P2(t) in black. 
The initial value of the first coordinate is larger than that of the second and 
third coordinates, and one can easily see from the graphs that our theoretical 
predictions are in full accordance with the properties represented by the fig- 
ure: the solutions are separated from zero by some positive numbers, coordinate 
Po{t) is always larger than Pi(t),P2{t), and we may suppose that the solution 
is approaching the discontinuous almost periodic solution of the system, as time 
increases. 
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Fig. 2 shows that there exists a time of stabihzation of a solution with initial 
value P{0) = (0.03, 0.5, 0.2) so that the first coordinate is smaller than the second 
and third ones at the initial moment t = 0. 




time t 



Figure 2: We can see that despite the "odd" relation between the coordinates 
at the initial moment, the normal state is eventually achieved, when coordinate 
Po{t) is larger than Pi{t) and P2{t)- The stabilization time of the solution is no 
more than 5 units. Again, we have the consistency of the theoretical result with 
the simulation observations. 
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4 Conclusion 



As follows from the preceding discussion, if system ([T]) satisfies certain conditions, 
then systemic arterial pressure and the blood pressure in compartments oscillate, 
remaining positive. The shape of the oscillations depends on the behavior of the 
moments when contraction of the left ventricle takes place. The oscillations 
are asymptotically stable. That is, they do not react significantly to external 
perturbations. 

From the results of our paper one can see that the strike type influence of 
the left ventricle contraction on the systemic arterial pressure is "softened" when 
it reaches the peripheral compartments through the connections of the system, 
so that the compartments' oscillations are continuous. On the other hand, this 
influence is sufficiently large to sustain the positiveness of all coordinates and 
keep the aortic pressure higher than in any other compartment. 

We believe that the proposed model can be developed further, using exper- 
imental results as well as the methods of the theories of functional differential 
equations, discrete equations, etc, to obtain additional features of the regular be- 
havior of blood pressure, as well as to investigate irregular (chaotic) processes in 
the system. The problem of the period-doubling and intermittency routes to the 
chaos can serve as an important theoretical tool in the studies of the arythmia 
processes. 
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